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Abstract 

A seismic model for the Sun calculated using the accurate helioseismic data 
predicts a lower neutrino flux as compared to the standard solar model 
(SSM). However, there persists a discrepancy between the predicted and mea- 
sured neutrino fluxes and it seems necessary to invoke neutrino oscillations 
to explain the measurements. In this work, we have performed a global, uni- 
fied oscillation analysis of the latest solar neutrino data (including the results 
of SNO charged current rate) using the seismic model fluxes as theoretical 
predictions. We determine the best-flt values of the neutrino oscillation pa- 
rameters and the Xmin both fe — inactive and Ve — t'steriie cascs and present 
the allowed parameter regions in the Am^ — tan^ 9 plane for fe — ^active tran- 
sition. The results are compared with those obtained using the latest SSM 
by Bahcall and his collaborators. 
PACS: 14.60.Pq, 12.15.Ff, 26.65.+t, 96.60.Ly 

I. INTRODUCTION 

Solar neutrino fluxes measured by all the experiments to date are significantly at variance 
with the expected theoretical predictions. The most recent confirmation of this has come 
from the heavy water (D2O) detector at Sudbury Neutrino Observatory (SNO) |Q which 
measures the solar neutrinos through the charged current (CC) reaction z/gd ppe~. 
SNO has also published their result of flux measured by the neutrino-electron scattering 
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reaction (ES) and reported a lower ®B flux as compared to theoretical predictions of the stan- 
dard solar model. This is in agreement with the flux measured by the SuperKamiokande 
(SK) detector through the same reaction 0|. Thus, SNO and SK conflrm the deflcit of solar 
neutrino fluxes reported flrst in the ^^Cl radiochemical experiment of Davis et al. [§] and 
subsequently by Kamiokande and the radiochemical ^^Ga experiments SAGE, GALLEX 
and GNO ||^. The theoretical predictions most widely used are from the standard solar 
model (SSM) developed and remarkably reflned over the last four decades by Bahcall and 
his collaborators In recent years, the observations of solar oscillations have provided 

an independent test of solar models. Inversions of accurately measured frequencies of solar 
oscillations have enabled a determination of the sound speed and density proflles inside the 
Sun P]. While the SSM matches these inverted proflles remarkably well, there is still a 
signiflcant discrepancy which is much larger than the errors in helioseismic inversions. It 
would, therefore, be desirable to check the results on neutrino oscillation solutions to the 
solar neutrino problem using a solar model which is consistent with the helioseismic data. 
Such models can be constructed using the inverted proflles of sound speed and density along 
with the equations of thermal equilibrium, provided the heavy element abundance proflle 
and the input physics like opacity, equation of state and nuclear reaction rates are assumed 
to be known [|T^,^. These models represent the present Sun and do not depend on evo- 
lutionary history of the Sun. Such seismic models can be used to calculate solar neutrino 
fluxes which turn out to be somewhat different from those obtained with the standard solar 
model. 

In this work, we consider the seismic model calculated using the technique described by 
Antia and Chitre but using updated helioseismic data. This model predicts a lower *B 



flux than the BP2000 SSM [Q. However, when all the experimental rates are taken together 
there is still inconsistency between theory and experiment |TT|. This inconsistency cannot 
be removed even if opacities and heavy element abundances are varied by arbitrarily large 
amount []T2[. Thus one needs to invoke neutrino oscillations to explain the observed fluxes 



of solar neutrinos. Seismic models have also provided some constraints on the pp reaction 



cross-section |Tl]JT^-|15| . It appears that pp reaction cross-section needs to be increased by 
about 4% over the value used by Bahcall et al. to obtain solar models that are consistent 
with seismic data. This increase in pp reaction rate in the seismic model with correct 
luminosity tends to decrease the predicted neutrino fluxes for all four experiments. 

In addition to the data on the total flux, SK also provides the data on the day-night 
recoil electron energy spectrum and the zenith angle distribution of events SNO 
has also reported the recoil electron energy spectrum for the CC events. Global oscillation 
analysis of the data on rates and spectrum has been carried out by different groups to 
put constraints on the oscillation parameters [lT8|-p3[] (pre-SNO) and p^^7| (post-SNO) 
These studies have used the neutrino fluxes from the standard solar model of Bahcall et. 
al. However, it would be interesting to see the implications of neutrino fluxes from the 

seismic model on the oscillation parameters. In this work, we perform a global and unifled 
oscillation analysis of the solar neutrino data on total fluxes from the SNO, CI, Ga and 
SuperKamiokande experiments together with the day-night spectrum data from SK using 
the seismically inferred neutrino fluxes and compare the results with those obtained using 
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the latest standard solar model (BP2000) of Bahcall et. al. []. We use the latest 1258-day 



SK data in our analysis [28 



The rest of the paper is organized as follows: in section 2 we present the basic features 
of the seismic model and describe the formalism for the unified analysis in Section 3. In 
section 4 we discuss the procedure for the analysis and the results and finally, summarize 
the conclusions in Section 5. 



II. THE SEISMIC MODEL 



We use the mean frequencies from the Michelson Doppler Imager (MDI) data from the 
first 360 days of its operation to calculate the sound speed and density profiles inside 



the Sun using a regularized least squares inversion technique [0. We adopt the inverted 
sound speed and density profiles, along with the heavy element abundance {Z) profile from 
the solar model of BTI] to calculate the seismic model using equations of thermal equilibrium 



rT| . We use the OPAL opacities |3^ and equation of state [Q, and nuclear reaction rates 
from [^], except for the proton-proton reaction rate, which is slightly adjusted to yield the 
correct observed solar luminosity With the latest physical and seismic input, the pp 



reaction cross-section needs to be increased to 4.17 x 10~^^ MeV barns. The helium 
abundance by mass at the base of the convection zone in this seismic model turns out to be 
0.250, which is in good agreement with the helioseismic estimate in the envelope . 



With the knowledge of temperature, density and composition profiles from the seismic 
model, neutrino fluxes can be calculated and the results are shown in Table |. Apart from 
inversion errors, other main sources of uncertainties in these calculated fluxes are the nuclear 
reaction cross-sections for the ^He-^He (5*3^3), ^He-^He (5*3^4), p-^^N (5*1^14), p-^Be {Sij) 
reactions, as well as the solar luminosity, the heavy element abundance, Z and opacities, n. 
To estimate the effect of these quantities on neutrino fluxes in the seismic model we calculate 
the logarithmic derivatives of neutrino fluxes with respect to each of these quantities {Xi) 
and these are also listed in Table ||, with the last row showing the estimated relative errors 
in these quantities. Apart from these we also include the uncertainty due to the electron 
capture rate of the process ^Be(e~, t'e)^Li and the astrophysical uncertainty in the 5*0 factor 
of the reaction ^^0(p, 7)^^F. These contributions are same as in |Q. The expected neutrino 
fluxes in various solar neutrino experiments can be calculated from Table |I| and these values 
are given in Table |I|, for comparison with observed fluxes 0,^,28| and those in the standard 



solar models P,[7|, |5B| , |57| . In Table |ITT| we present the contributions of the various neutrino 
sources to the CI and Ga experiments according to the seismic model and BP2000. Table 
|V] summarises the ratios of the experimental rates to the theoretical predictions for the CI, 
Ga, SK and SNO experiments for both BP2000 and seismic model. For SNO we display 
only the CC rate. We also show the compositions of major flux components. 

The neutrino fluxes in our seismic model are somewhat different from those in the seismic 



model of Watanabe and Shibahashi [38]. The main difference arises because they have used 



We have not incorporated the SNO ES data as well as the SNO CC spectrum data as they still 
have large errors. 
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only the sound speed from primary inversions to calculate the seismic model, while we have 
used both sound speed and density profiles from primary inversions. Thus the density profile 
in seismic model of Watanabe and Shibahashi does not, in fact, match the inverted density 
profile. Further, since the density profile in their model is not constrained to seismic profile 
they are able to construct seismic model with correct luminosity using the standard pp 
nuclear reaction rate, which tends to give larger neutrino fluxes. If they were to assume 
a slightly larger cross-section for pp reaction, the density profile in their model would be 
in better agreement with the inverted profile. Apart from this, the adopted heavy element 
abundance, Z-profile is also different in their work. While we have taken the Z-profile from 



model of Richard et al. |31|, which includes mixing just below the base of convection zone. 



as implied by helioseismic data Watanabe and Shibahashi use a homogeneous Z 



profile. Similarly, the neutrino fluxes in our seismic model are somewhat smaller than those 
in SSM [§]. The main reasons for this reduction are again the increase in S\\ and neglect of 
mixing just beneath the convection zone in SSM. If these had been incorporated, then the 
SSM fluxes should also come close to the corresponding seismic model fluxes. 

III. FORMALISM FOR UNIFIED OSCILLATION ANALYSIS 

The general expression for the probability amplitude of survival for an electron neutrino 
produced in the deep interior of the Sun, for two neutrino flavors, is given by |^ 

A — A® /ivac 4® I 40 /ivac /I® 

where Afi^{k = 1, 2) gives the probability amplitude of u^. — > z/^ transition at the solar surface, 
AJ.^ is the survival amplitude from the solar surface to the surface of the Earth and A®^ 
denotes the Uk — ^ transition amplitudes inside the Earth. We can express 

A% = a%e-^'^^ (2) 

where (pf is the phase picked up by the neutrinos on their way from the production point 
in the central regions to the surface of the Sun and 

af,' = l + il-Pj)cos2e^ (3) 

9m being the mixing angle at the production point of the neutrino, given by 

Am^sin2^ 

tan 2^^ = ^ . (4) 

Am^ cos 29 -2 V2GFneE ^ ' 

Here rif, is the ambient electron density, E the neutrino energy, and Am^ (= — m^) is 
the mass squared difference in vacuum. We denote by Pj the non-adiabatic level jumping 
probability between the two mass eigenstates which for an exponential density profile^ can 
be expressed as El 



^Note that for the actual calculation of the probabilities we have used the numerical density profile 



given in |8[ for SSM and the model in |11| with updated helioseismic data for seismic model. 
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exp(— 7 sin^ 6) — exp(— 7) 
1 - exp(-7) 



(5) 



Am? 



7 



TT- 



E 



dr 



For sterile neutrinos, the in Eq. has to be replaced by Ue 
neutron number density. The survival amplitude A'j^ is given by 



(6) 



|n„, where n„ is the 



A 



kk 



-iEk{L-RQ) 



(7) 



where Ek is the energy of the state Uk, L is the distance between the center of the Sun and 
Earth and Rq is the solar radius. For a two-slab model of the Earth — a mantle and core 
with constant densities of 4.5 and 11.5 gm cm~^ respectively, the expression for Af^ can be 
written as (assuming the flavor states to be continuous across the boundaries) 



a: 



2e 



E 

a,l3,(T 



ak ai 



(8) 



where {i,j,k) denotes mass eigenstates and (a,/?, cr) denotes flavor eigenstates, U, and 

are the mixing matrices in vacuum, in the mantle and the core respectively and if) 
and ip^ are the corresponding phases picked up by the neutrinos as they travel through the 
mantle and the core of the Earth. The Ve survival probability is given by 



P - 14 P 



„02| 4® |2 , ^02| 40 |2 
"el l^lel ^ "e2 l^2el 



i{02,0-<^i,©)l 



Identifying = a^i and P© = |v4® p Eq. (H) can be expressed as [|3],|4[|0 



Pee = P0Pe + (l-P0)(l-P0 



+2v/P0(l-P0)Pe(l-P0)cose 



(9) 

(10) 
(11) 



where we have combined all the phases involved in the Sun, vacuum and inside Earth in ^. 
This is the most general expression for survival probability for the unified analysis of solar 
neutrino data. Depending on the value of Am? / E one recovers the well known Mikheyev- 
Smirnov-Wolfenstein (MSW) and vacuum oscillation (VO) limits: 

• In the regime Am'^/E ~ 5 X 10 eV^/MeV matter effects inside the Sun suppress 
flavor transitions and 6m ~ 7c/2. Therefore, from (^), we obtain Pq ^ Pj ^ cos^ 6 as 
the propagation of neutrinos is extremely non-adiabatic and likewise, P© = cos^ 6 to 
give 



1 - sm'^29sm^{Am^ (L - Rq)/4:E) 



(12) 
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For ^m? / E ~ 10~^ eV^/MeV, the total oscillation phase becomes very large and the 
cos^ term in Eq. (|TID averages out to zero. One then recovers the usual MSW survival 
probability 

pMsw _pD^ (2P^-l)(sin^g-|A|p) 

The day-time probability is given by being 

= ^ + - Pj) cos 29 cos 26^ (14) 

In between the pure vacuum oscillation regime where the matter effects can be safely 
neglected, and the pure MSW zone where the coherence effects due to the phase ^ 
can be conveniently disregarded, is a region where both effects can contribute. For 
5 X 10"^° eVVMeV ~ Arn^ / E ~ 10"^ eVVMeV, both matter effects inside the Sun 
and coherent oscillation effects in the vacuum become important. This is the quasi 
vacuum oscillation (QVO) regime In this region, Pq ^ Pj and P^ = cos^ 6' and 



the survival probability is given by [^,^ 



Pee = Pj cos^ ^ + (1 - Pj) sin^ 9 + sin^ 20^Pj{l-Pj) cos^ (15) 
we calculate Pj in this region using the prescription given in M7 . 



IV. RESULTS AND DISCUSSIONS 

In this section we present our results of global -analysis of the data on total rates and 
the day-night spectrum observed in SK. We have done two sets of calculations taking the 
theoretical predictions from the seismic model as well as from BP2000 SSM. Our principal 
objective is to compare the two sets of results. 

Different approaches have been adopted regarding the treatment of the SK data on total 
rates and the recoil electron energy spectrum, in the global analysis. A critical comparison of 
the various methods used and the dependence of the final results on the method followed is 
discussed in a lucid and extensive manner in [^. Below we summarise the salient features. 



(i) One approach is to to include both the SK rate and the spectrum in the global 
analysis [ PT| , P^ . 

(ii) The total rate measured in SK is not independent of the spectrum and hence the 
above method of including both may lead to a possible overcounting of events. To 



avoid this the total SK rate is excluded from the global analysis in p2| , p5 



(iii) Another approach adopted to avoid the overcounting is to include the total SK 
rate in the global analysis but to adopt a free normalization factor for the spec- 



trum [|19| , p2| - |2^j26|] . Thus the spectrum gives information on only the shape of the 
survival probability. 
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In case of methods (i) and (ii), if the flux normahsation is varied as a free parameter 



in both rates and spectrum, no error correlation between these needs to be taken |2^. On 
the other hand, if the flux normalisation is kept fixed at the SSM value, then one should 
incorporate the correlation between the error in the the ^B flux measured in the total rates 



and the SK spectrum ||2l|j25[1 . For method (iii), since the normalisation of the spectrum is 
varied as a free parameter, the due to rates and spectrum can be summed independently. 
For the purpose of this paper we adopt method (iii), so that our total is given by 

Xtotal Xrates ~^ Xspectrum 

(16) 

where 

xLes = E (Pf - {<y^) {Ff^ - ■ (17) 

Here = T//Tf^^, where ^ designates th (for the theoretical prediction with oscillations) 
or exp (for the experimental value) and Tj stand for the quantities total rates from different 
experiments. We have used the weighted average of the three Ga experimental rates and 
thus we have 4 data points for the total rates. The error matrix cxjj contains the experimental 
errors, theoretical errors and their correlations. The correlation matrix for the total rates 



is constructed as in The logarithmic derivatives for the seismic model needed for the 
calculation of error correlation matrix are given in Table |. For the SSM we take the latest 
values from [||. The spectrum chi-square is defined as, 

Xspectmm = E (^sp-^*"" " ^T'''') {<^i'/) {^spFj^^ " i^""^) • (18) 

Here = Sf / Sf^^ where ^ designates th (for the theoretical prediction with oscillations) or 
exp (for the experimental value) and Si stands for the electron energy spectrum for different 
energy bins. We have used 1258 day SK data on the electron energy spectrum at day and 
night which includes the energy bin from 5.0-5.5 MeV also and we have 38 data points for 
the spectrum. Xgp is the normalisation factor for the spectrum which is floated as a free 
parameter in the global analysis in order to filter out the information on the total flux from 
the spectrum data. Thus for the global analysis we have 41 degrees of freedom (DOF) for the 
no oscillation scenario. For the error matrix in the spectrum data we include the statistical 
error, correlated and uncorrelated systematic errors and the error due to the calculation of 
the spectrum |T§. The no oscillation xLn 41 DOF for BP2000 SSM is 89.27 while for 
the seismic model it is 94.17. 

Next, we perform the chi-square analysis assuming neutrino oscillations to operate. For 
the oscillation analysis, there are two parameters - Am^ and tan^ 6 and thus the number of 
DOF is 39. In Table |V| we present the best-fit values of parameters, xLin goodness 
of fit (GOF) of the solutions for both BP2000 and seismic model for z/g — t'active transitions. 
There are five allowed solution in both models - Large Mixing Angle (LMA), Vacuum 
Oscillation (VO), Low Am^ - Quasivacuum Oscillation (LOW-QVO), Just So^ [^] and 
Small Mixing Angle (SMA) - in order of decreasing GOF. The GOF in these regions are 
more or less similar in both models. The best-fit for both models comes in the LMA region. 
In order to understand the results of Table 0, we write Eq. (|13D as 
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Pee = P''-\{2P''-l)fre, (19) 

with e = cos 26', /reg = |^®eP — sin^6' and P^ defined above in Eq. (p^. In Fig. 1 we plot 
the regeneration factor /reg and the actual Earth regeneration Re = Pee — vs. energy at 
the SK latitude for the best-fit values of parameters in the SMA, LMA and LOW regions. 
Since the latitude of the other detectors are not very different we do not expect /reg and Re 
to be very different for these. Since /reg is always positive, the possibility of regeneration 
inside the Earth depends on Pj For Pj < | there is u^, regeneration inside Earth while for 
Pj > \ more z/g's are flavor converted. 

• For the SMA region e ~ 1 and from Fig. 1 we observe that /reg is very small excepting 
for two peaks at E ^ 6 MeV and E ^ 15 MeV corresponding to strong enhancement of 



the earth regeneration effect for the neutrinos passing through the core p9| , p2| . Hence 

pSMA ^ pD (20) 



In this region for low energy {pp) neutrinos, resonance is not encountered (reso- 
nance density ^ maximum solar density) and hence Pj ~ and cos 29m ~ 1 giving 
Pee^^ ~ 1. For intermediate energy (J Be) neutrinos cos2^m ~ —1 (resonance density 
■C production density) and P^^^ Pj fa for these energies. For high energy (^P) 
neutrinos also, cos2^m ^ — 1 and P^^^ ~ Pj, with Pj rising with energy. This energy 
dependence of the SMA survival probability gives a GOF of 9.22% for BP2000 and 
9.21% for the seismic model for a simultaneous description of all the four observed 
rates given in Table |IV| and the SK spectrum. 

• For the Am^ of the LMA solution in Table 0, the motion of the neutrino in the 
solar matter is adiabatic for almost all neutrino energies and Pj ~ 0. For low energy 
neutrinos the matter effects are weak both inside the Sun and in Earth giving /reg ~ 
and cos 29^ ~ e so that for Ga energies 

Pi;^^^i(l + e2) (21) 



At SK and SNO energies matter effects result in cos26'm ^ — 1 while /reg is small but 
non-zero (~ 0.03 at 10 MeV as seen from fig. 1) giving 

n"e^^^^(l-e) + /reg 

= Sin^ 9 + /reg (22) 



With the values of e from Table and /reg given in fig. 1 eqs. (|2ll) and (^2|) approx- 
imately reproduce the rates of Table Since the probability p^ ) is approximately 
energy independent it can account for the fiat recoil electron energy spectrum. Since 
the seismic model predicts a higher value for observed to predicted ratio for the SK, 
SNO and CI rates, the best-fit value of sin^ 9 or tan^ 9 obtained for seismic model are 
larger (cf. Eq. (^)). 
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• In the LOW region cos 26m ~ — 1 for all neutrino energies and Pj ~ (except for very 
high energy neutrinos) and thus 

n™ = ^(l-e)+Aeg (23) 

As is seen from Fig. 1 /j-eg is small for high energy neutrinos and large for low energy 
neutrinos. For the best-fit case e = 0.2 and /reg ~ 0.2 for Ga energies and ~ 0.025 for 
SK energies, which can just about reconcile the Ga and SK rates. But it provides a 
very good description of the flat SK spectrum. The best-fit mixing angle is larger for 
the seismic case as in LMA. 

• At the best-fit value of Am^ of Table for VO solutions the energy smearing over the 
bins washes out the energy variation due to the oscillations and the fiat recoil electron 
spectrum observed at SK can be explained. 

• For the Am^ in the Just So^ region one gets a very small survival probability for the 
^Be neutrinos while for the neutrinos the survival probability is close to 1.0 [|52| . 
Therefore it cannot explain the total rates data but since it gives a fiat probability 
for the *B neutrinos the spectrum shape can be accounted for and the global analysis 
gives a GOF of 8.1% in BP2000 and 12.49% in seismic. Since the ratios of observed 
rates to predicted rates are higher for seismic Just So^ gives a lower contribution to 
Xrates a better GOF in seismic. 

In Fig. 2 we plot the 90% (x' < xLn + 4.61), 95% {x' < xLn + 5-99), 99% {x' < 
Xrnin + 9.21) and 99.73% C.L.(x^ < xLn + 11-83) allowed areas in the Am^ - tan^^ plane 
for both BP2000 and the seismic model. We plot the allowed regions with respect to the 
global minimum. Both the models do not admit any allowed area in the SMA region. For 
other regions the allowed areas in SSM and seismic model are roughly similar with the 
following differences observed in the seismic case: 

• more area seems to be allowed in the LMA and LOW-QVO region 

• the LMA region extends to higher Am^ 

• higher values of mixing angles are allowed for LMA and LOW-QVO regions, specifi- 
cally, with LMA extending into the dark zone {6 > vr/4). 

For the Ga experiment, the seismic model predicts a higher pp and lower flux as 
compared to BP2000. The net effect is that the increase of the flux ratio is smaller in Ga 
than in SNO, SK and CI and the observed to predicted rate in Ga is closer to that in CI, 
SK and SNO for the seismic model. Thus the energy dependence between the low energy 
pp and high energy ^B rates seen for BP2000 is reduced in the seismic model and the data 
can now be better accounted by an energy independent scenario. In fact, since the energy 
distortion in the observed rate is less in the seismic case, most parts of the bands in the 
parameter space with energy dependence of < 10% are permitted^. This accounts for the 



^ Figure 1 of reference |^ . 
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extended allowed areas appearing in Fig. 2 for the seismic model. These extended areas 
with weak energy dependence include both high Am^ zones as well as high mixing angle^. 
In fact, the allowed area expands well into the dark zone for the LMA solution. For these 
values of the mixing angles the predicted flux ratio for the SK and SNO is more than that 
for Ga, against the energy trend of the data. However, for the seismic model these zones 
are still allowed at 99.73% CL. owing to the proximity of the Ga and SK-SNO rates. 

In Fig. 3a and Fig. 3b we show the experimental rates with their 3cr errors, together 
with the 3cr range of predicted rates for the LMA (Fig. 3a) and LOW-QVO-VO (Fig. 3b) 
solutions, in the plane of any two experiments taken together. From the figures it is clear 
that the LMA (LOW-QVO) region can better account for the SSM(seismic model) predicted 
rates. 

In Table |V| we give the xLin ^^e best-fit points for Ue transition to sterile neutrino. 
Since the flux measured by the charged current reaction at SNO is significantly lower 
than the observed SK flux, all the sterile solutions appear to be disfavored with more than 
90% probability, except for the VO solution for BP2000 and both VO and Just So^ for the 
seismic model. The VO solution produces better fits than the MSW solutions as it gives a 
lower contribution to the Xspectrum (cf.eq. (|l^). As in the z/g — inactive case. Just So^ gives a 
better fit for seismic. 



V. CONCLUSIONS 

The measured solar neutrino fluxes have been consistently lower than the theoretical 
predictions from SSM. We have constructed a seismic model of the sun consistent with 
the helioseismic data that predicts fluxes lower than that predicted by BP2000, but the 
corresponding pp flux turns out to be slightly higher in the seismic model. We examine 
how the use of the seismic model fluxes changes quality of the fits in the MSW and the 
vacuum oscillation region as compared to BP2000. For the statistical analysis of the data 
we use a minimization technique where we vary the normalization of the spectrum as a 
free parameter and thus avoid the over-counting of the SK observed flux and consider only 
the shape information from the SK spectral data. We find that the inclusion of theoretical 
uncertainties and their correlations self consistently in both models result in fairly similar 
GOF for the oscillation solutions in both models. However, we note that the use of seismic 
fluxes does modify the allowed areas in the parameter space. The increased proximity of the 
Ga and SK rates reduces the energy distortion of the observed fluxes and allows the regions 
of parameter space with weak energy dependence. 
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TABLES 



TABLE I. Neutrino fluxes in seismie model 



Neutrino 


Flux, Fy 

(cm s j 


'S'3,3 


Lo^ 

'S'3,4 


jarithmic derivatives g'^'^^ 

'Si, 14 'Si, 7 Lq 


Z 


K 




f6 12 ± fl Oil X 10^° 


0.03 


—0.07 


-0.02 


0.00 


0.96 


—0.06 


—0.13 


pep 


(1.45 ±0.01) X 10^ 


0.03 


-0.07 


-0.02 


0.00 


0.76 


-0.17 


-0.31 


hep 


(2.12 ±0.01) X 10^ 


-0.41 


-0.07 


-0.01 


0.00 


0.09 


-0.21 


-0.41 


^Be 


(4.54 ±0.12) X 10^ 


-0.44 


0.95 


0.00 


0.00 


2.06 


0.68 


1.45 


8b 


(4.16 ± 0.26) X 10^ 


-0.49 


1.04 


0.01 


1.00 


3.85 


1.65 


3.40 


13n 


(5.02 ± 0.25) X 10^ 


-0.08 


0.14 


0.85 


0.00 


2.39 


1.24 


2.52 


150 


(4.14 ±0.26) X 10^ 


-0.09 


0.17 


1.00 


0.00 


2.93 


1.52 


3.07 


17f 


(4.77 ± 0.32) X 10^ 


-0.10 


0.18 


0.01 


0.00 


3.07 


1.59 


3.21 


Rel. error 




0.060 


0.094 


0.143 


0.11 


.004 


0.033 


0.02 



TABLE II. Predicted neutrino fluxes in various solar neutrino experiments 



Experiment Homestake SK Gallex, SAGE, GNO SNO CC 

(37C1) (H2O) C'Ga) (D2O) 

(SNU) (lO^cm-^s-i) (SNU) (lO^cm-^s'i) 

Measured Flux 2.56 ± 0.23 2.32 ± 0.08 74.7 ±5.0 1.75 ±0.14 

SSM (BPOO) 7.6li:? 5.05lJ;^} 128.0l? 5.05^^;^^ 

SSM (BP95) 9.3lj:4 6.621?;?^ 1371? 6.621^3 

SSM (BP98) 7.7li:g 5.151^;^^ 129^^ 5.151^;^^ 

SSM (Bru98) 7.2 4.8 127 4.8 

SSM (Bru99) 6.7 4.7 125 4.7 

Seismic model 6.46 ± 0.99 4.16 ±0.76 124.9 ±6.5 4.16 ±0.76 
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TABLE III. 

The seismic model predictions for the solar neutrinos fluxes and neutrino capture rates in 
the CI and Ga detectors. The expected ®B flux in SK and SNO is as given in Table |1[ 





bcloilllC 




RP9nnn 




oU Ui 




71 Ha 




\sxcX, 








I^OIN KJ j 




pp 


0.00 


71.58 


0.00 


69.7 


pep 


n OQ 
U.zo 


z.yt) 


U.zz 


z.o 


hep 


U.U4 


U.U ( 


n n/i 

U.U4 


n 1 

U.l 


^Be 


1.09 


32.57 


1.15 


34.2 


8B 


4.74 


9.98 


5.76 


12.1 


13n 


0.09 


3.03 


0.09 


3.4 


150 


0.27 


4.72 


0.33 


5.5 


17p 


0.00 


0.05 


0.00 


0.1 


Total 


6.46 


124.94 


7.6 


128.0 



TABLE IV. The observed rates relative to the theoretical predictions for Ga, CI, SK and SNO 
experiments along with their compositions(comp.) for both seismic model and SSM. The Ga rate 
corresponds to the combined SAGE and GALLEX+GNO data. 



Model 


Experiment 


Ga 




CI 


SK 


SNO (CC) 




Rate 


0.584 ±0.039 




0.337 ±0.030 


0.459 ±0.017 


0.347 ±0.027 


BP2000 
















Comp. 


pp (55%), ^Be (25%), ^ 


*B (10%) 


^Be (15%), 8b (75%) 


SB (100%) 


SB (100%) 




Rate 


0.598 ±0.040 




0.396 ±0.035 


0.557 ±0.021 


0.421 ±0.033 


SEISMIC 
















Comp. 


pp (57%), ^Be (26%), 


«B (8%) 


^Be (17%), (73%) 


SB (100%) 


SB (100%) 
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TABLE V. The best-fit values of the parameters, Xmin' goodness of fit from the global 

analysis of rates and the 1258 day SK day-night spectrum data for MSW analysis involving two 
active neutrino flavors. 





Nature of 




tan2 6l 


Xmin 


Goodness 




Solution 


in eV^ 






of fit 




SMA 


5.28 X 10"^ 


3.75 X 10-^ 


51.14 


9.22% 




LMA 


A Tn 1 n — 5 
4. lU X Lu 


U.60 


66AZ 


i 2.18/0 


BP2000 


LOW-QVO 


1.76 X 10"^ 


0.67 


39.00 


46.99% 




VO 


4.64 X 10-1° 


0.57 


38.28 


50.25% 




Just So^ 


5.37 X 10-12 


0.77 


51.90 


8.10% 




SMA 


4.66 X IQ-s 


5.10 X 10-^ 


51.15 


9.21% 




LMA 


5.11 X 10-^ 


0.44 


35.62 


62.48% 


SEISMIC 


LOW-QVO 


1.76 X 10-"^ 


0.71 


38.22 


50.53% 




VO 


4.65 X 10-1° 


0.47 


36.23 


59.69% 




Just So^ 


5.37 X 10-12 


1.00 


49.30 


12.49% 


TABLE VI. 


The The best-fit values of the parameters, Xmin' ■ 


goodness of fit from the 


global analysis of rates and the 1258 day SK day-nij 


^ht spectrum data for two-generation 


^sterile 


MSW analysis. 














Nature of 




tan2 6l 


A mm 


Goodness 




Solution 


in eV^ 






of fit 




SMA 


5.59 X lO-'' 


2.83 X 10-^ 


54.21 


5.35% 




LMA 


6.13 X 10-5 


0.50 


52.93 


6.75% 


BP2000 


LOW-QVO 


2.93 X 10-^ 


1.00 


53.18 


6.46% 




VO 


4.67 X 10-1° 


0.37 


46.28 


19.70% 




Just So^ 


5.37 X 10-12 


0.77 


52.09 


7.83% 




SMA 


3.81 X 10-6 


3.67 X 10-^ 


58.18 


2.47% 




LMA 


6.04 X 10-5 


0.62 


54.97 


4.64% 


SEISMIC 


LOW-QVO 


3.20 X 10-^ 


1.00 


53.26 


6.36% 




VO 


4.68 X 10-1° 


0.37 


44.82 


24.09% 




Just So^ 


5.37 X 10-12 


0.98 


49.51 


12.07% 
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Fig. 1: The regeneration factor /reg = (|A2ep — sin^ 9) and the net Earth regeneration 
Re{= Pee — P^) as a function of energy for typical values of the parameters in the 
SMA, LMA and LOW-QVO regions. 
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tan e tan 

Fig. 2: The 90%, 95%, 99% and 99.73% C.L. allowed areas from the global analysis 
of the total rates from CI, Ga, SK and SNO detectors and the 1258 day SK recoil 
electron spectrum at day and night, assuming MSW conversions to active neutrinos, 
using theoretical predictions from BP2000 and seismic model. 
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Fig. 3a: The experimental rates (in SNU for the CI and Ga experiments and in 

units of 10^ cm~^s~^ for SK and SNO) with Scr errors (shown by ehipses) and the 
99.73% C.L. range of predicted rates for LMA solution, in the plane of any two of the 
experiments for BP2000 (dotted line) and for seismic model (solid line). 
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Fig. 3b: same as in fig. 3a but for LOW-QVO-VO solution. For this solution there 
are multiple contours (cf. Fig. 2), resulting in a complicated pattern. 
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